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Abstract 



We consider quantum systems consisting of a linear chain of n harmonic oscillators coupled 
by a nearest neighbour interaction of the form —qrcjr+i {cjr refers to the position of the rth 
, oscillator). In principle, such systems are always numerically solvable and involve the eigenvalues 

of the interaction matrix. In this paper, we investigate when such a system is analytically 
(~| . solvable, i.e. when the eigenvalues and eigenvectors of the interaction matrix have analytically 

closed expressions. This is the case when the interaction matrix coincides with the Jacobi matrix 
of a system of discrete orthogonal polynomials. Our study of possible systems leads to three 
new analytically solvable Hamiltonians: with a Krawtchouk interaction, a Hahn interaction or 
a (7-Krawtchouk interaction. For each of these cases, we give the spectrum of the Hamiltonian 
(in analytic form) and discuss some typical properties of the spectra. 

> 

00 . 

^ ' 1 Introduction 
(N , 

I In classical mechanics, one-dimensional systems (or lattices) consisting of mass points with some 
■ nearest neighbour interaction have a long history. A typical system is a lattice of n particles with 
On . masses mi, 1712, ■ ■ ■ , rUn, and a harmonic coupling with spring constants ki, K2, ■ ■ ■ , Hn-i leading to 
' the Hamiltonian 

n 2 n— 1 

•O: ^(^'>-) = E^ + Ey(-.-W- (1) 

. r=l ^ r=l 

S ' • • 

Such classical systems (or variations, with an infinite number of mass points, or with various 
boundary conditions) were already considered by Schrodinger [1]. The equations of motion of 
such a system can be solved by (numerically) diagonalizing the interaction matrix (the eigenvalues 
of which yield the normal modes of the system). Alternatively, the system can be solved using 
orthogonal polynomials whose recurrence relations are derived from the equations of motion [2-4]. 
In that case, the normal modes are obtained from the zeros of the nth orthogonal polynomial. 

One-dimensional systems with a different type of nearest neighbour interaction received a lot of 
attention, especially those that are still exactly solvable. Among the most famous, we mention the 
Toda system [5] and the Calogero-Sutherland-Moser models [6-8]. In this context, the emphasis 
was shifted from physics to mathematical aspects such as integrability and the underlying algebraic 
structures. 

Also the quantum versions of many of these systems or models were investigated from various 
points of view during the last decades. Quantum Calogero-Moser systems for any root system were 
studied by Olshanetsky and Perelomov [9]; for a review, see [10]. In such quantum systems, the 
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emphasis - from the physics point of view - is on a construction of ground wave states, formulae 
for the excitation spectrum, a description of stationary states, etc. Several other quantum systems 
with a nearest neighbour interaction closely related to Calogero-Sutherland-Moser models were 
explored, see [11-14] to cite a few. 

In the present paper we consider yet another quantum system given in equation ([5]), consist- 
ing of a one-dimensional chain of particles with a certain nearest neighbour interaction, which is 
quadratic in the position operators. Our emphasis is on the investigation of analytical solvability 
of the quantum system (i.e. on obtaining analytically closed expressions of the spectrum of the 
Hamiltonian and of its eigenstates, the stationary states of the system). In such a context, the 
physical significance of the interaction introduced here is less clear: it can, in a sense, be considered 
as a deviation from a vibrating quantum system. 

To introduce our system, let us first consider one of the most common quantum systems, 
consisting of a chain of harmonic oscillators coupled by some nearest neighbour interaction [17- 
20]. In this popular model the particles are described as identical harmonic oscillators which are 
moreover coupled by springs obeying Hooke's law. Then the Hamiltonian of the system is given 
by: 

r=l ■r=0 

where % = Qn+i = (fixed wall boundary conditions). In other words, the quantum system consists 
of a string or chain of n identical harmonic oscillators, each having the same mass m and natural 
frequency oj. The position and momentum operators for the rth oscillator are given by and pj.] 
more precisely qr measures the displacement of the rth mass point with respect to its equilibrium 
position. The last term in ([5]) represents the nearest neighbour coupling by means of "springs", 
with a coupling strength c (c > 0). For c = one is simply dealing with a set of identical uncoupled 
harmonic oscillators. 

We shall consider ([21), and the other quantum system being studied here, in the case of canonical 
quantization, i.e. the qr and pr are self-adjoint operators (gj = qr and p\ = pr) satisfying the 
commutation relations 

[qr,qs]=0, [pr,ps]=0, [qr,Ps] = ih6rs (r, s = 1 , . . . , n) . (3) 

It is well known that the Hamiltonian ^ is completely solvable (see also next section). In fact, it is 
analytically solvable in the sense that one has an analytically closed expression for the eigenvalues 
of 

Note that ([2]) can be rewritten in the following form (with 0)^ = 0;^ + 2c) 

n „2 ~2 n— 1 

= E (£ + ^^Ir) -CmJ2 4rqr+l. (4) 
r=l r=l 

The fact that the nearest neighbour interaction is the same everywhere in the chain, thus indepen- 
dent of the position r, implies that all coefficients of ^r^r+i are the same. 

In the present paper, we shall study deviations of dH), where the nearest neighbour interaction 
is not constant in the chain, but depends on r. The Hamiltonian of such a more general system is 
assumed to be of the following form: 

n „o 9 n—l 

sr^ ( Vr cm v-^ ^ „ , , 

^ = E + - ^ E ^r^r^r^- (5) 

r=l r=l 

The physical interpretation of such a system is not obvious. First of all, in order to have a physical 
meaning, the interaction matrix related to ([5]) should still be positive definite (see next section). 



2 



Second, the interaction is in general no longer harmonic. It could be seen as a chain still consisting 
of identical oscillators but with a quadratic nearest neighbour interaction of the form — g^^r+i- This 
interaction is not homogeneous in the chain but depends on the location r in the linear system. 
For a more general context in which quantum systems of the form ([5]) appear as a special case, 
see the notion of "harmonic systems on general lattices" in [15,16] in the study of entanglement in 
many-body systems. 

The purpose of this paper is to study the spectrum of Hamiltonians of the form Such 
Hamiltonians are always numerically solvable, i.e. the complete spectrum can be described using 
the eigenvalues of the interaction matrix (see next section). The contribution of this paper is 
to investigate those cases for which H is analytically solvable (i.e. when the eigenvalues of the 
interaction matrix are known in analytically closed form). Some examples of Hamiltonians for 
which this is the case are: 

n ^9 9 n— 1 

r=l r=l 



rrii) \^ f Pr , "^^^ ■^2\ cm ^ \/(n - r) (n + r + 1) 



\2m 2 ^"^J 2 ^ 2 

■r=l r=l 



QrQr+1, 



n ^2 2 n— 1 

r=l r=l 

where q is some positive parameter. 

In the following section, we shall consider ^ again, and describe a method to solve this Hamil- 
tonian. This method can best be formulated in terms of an interaction matrix M [11,21], and 
we shall determine the conditions for M to be analytically solvable. In particular, systems of the 
form ([5]) correspond to tridiagonal interaction matrices with a constant diagonal. As we shall see, 
this leads us to the area of discrete orthogonal polynomials. We have investigated the families 
of discrete orthogonal polynomials that lead to solutions. These are described in the following 
sections: an interaction based upon Krawtchouk polynomials, on Hahn polynomials, or on dual 
g-Krawtchouk polynomials. For each of these cases, we give the corresponding Hamiltonian and 
its solution. In section [6] we also briefly describe some interesting features of the energy levels of 
these Hamiltonians. 

2 General method 

The Hamiltonian ([2]) can be written in matrix form as follows: 



( Pi 



+ f ( • • • qi 



( 

\ 9n 



(6) 



V Pr. 

where Ah is a symmetric tridiagonal matrix of the form 

Ah = uj'^I + cMh, (7) 
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with I the n X n identity matrix and 



/ 2 -1 

-1 2 -1 
0-12 



V 



\ 





2 -1 

-1 2 / 



(8) 



The fact that the interaction is only between nearest neighbours is reflected by the tridiagonal form 
of Mfi (i.e. nonzero entries only on the diagonal, the subdiagonal and the super diagonal). We shall 
refer to Mh as the "interaction matrix" [11,21]. 

Let us for a while consider quantum systems described by a Hamiltonian with a more general 
interaction matrix M: 



H 



1 
2m 



Pn 



( Pi 



\ Pn 



qt )iu;''I + cM) 



(9) 



In ([9]) , M is a real and symmetric matrix. In order to be physically meaning ful, uj'^I + cM should be 
a positive definite matrix [15,16]. A general method to deal with such Hamiltonians was described 
in [15, section 2.1]. Since M is real and symmetric, the spectral theorem [22] implies 



where 



M = UDU^ 



D = diag(Ai, A2, . . . , A„), 
uu^ = U^U = I. 



(10) 



(11) 

(12) 



The entries of the diagonal matrix D are the (real) eigenvalues Aj of M, in some order, and the 
columns of the real orthogonal matrix U are eigenvectors of M (in the same order); U'^ stands for 
the transpose of U. 

Introducing new operators (the so-called normal coordinates and momenta) as follows: 



the Hamiltonian 



H 



1 

2m 



f Qi 

\ Qn 

reads 

••• 



( 4i 




U 



T 



( Pi 



\ Pn 



(13) 



/ A 



m , 



qI (uj^I + cD) 



\ Pn 



( Qi\ 



\Qn ) 



2m 



(14) 



3=1 j=i 

By the transformation ()13p . the new operators also satisfy the canonical commutation relations: 
[Qj,Qk] = 0, [Pj,Pk]=0, [Qj,Pk]=ih6jk {j,k = l,...,n). (15) 



In (|14p . the values of iJ^ + cAj are all positive since the interaction matrix uP'l + cM is assumed to 
be positive definite. So one can introduce 



'a;2+cAj. (16) 
Then we can write 

* = 2^E^/ + ?E'^J<5?. (17) 

i=i i=i 

This expression is just like the Hamiltonian of an n-dimensional non-isotropic oscillator, so we can 
use the commonly known method for its solution [23,24]. Introducing boson operators 



these satisfy 

[a-,o^] = [a+,a+] = 0, [aj ,a'^] = 6jk, {j,k = 1, . . . ,n) (19) 
and H can be written as 

Furthermore, 

[H,af]=±nwjaf (i = l,...,n). (21) 

So if we assume that there is a lowest ^f-eigenvalue (lowest energy) , say for the state 1 0) , then we 
have the usual n-boson Fock space in which the action of H is diagonal. The vacuum vector |0) 
satisfies 

(0|0) = 1, a-|0) = 0; (22) 
the other (orthogonal and normalized) basis vectors are then defined by 



. . . , = |o), (kj = 0, 1, . . .). (23) 



The spectrum of H is now determined by 

" 1 

H\ki,...,kn) = '^hujj{kj + -) \ki,...,kn). (24) 
i=i 

This analysis is well known, and it seems to indicate that a Hamiltonian of the form ([9]) with 
a general interaction matrix M is exactly solvable as a quantum system. Note, however, that 
the solution we have described involves also a numerical process, namely the determination of 
the eigenvalues and eigenvectors of M in (jlip and (jl2p . We shall say that the Hamiltonian H is 
analytically solvable if we have an analytically closed expression for the eigenvalues and eigenvectors 
of M, for arbitrary n. 

One example of an analytically solvable Hamiltonian is ([2]), with interaction matrix AI^ given 
by dS]). In this case, the decomposition (fTO|) is determined by [17,20] 

C/=^(si„(^)) (25) 
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and the eigenvalues in pT]) by 

A,- = 2 - 2cos(^^). (26) 

So we have 

J] = u? ^2c- 2ccos(^^) = + 4csin^( J . (27) 

n + 1 2(n + 1) 

Apparently, there are not so many examples of analytically solvable Hamiltonians of the type Q 
known in the literature. One paper dealing with this problem (and closely related ones) is [21]. In 
that paper, some examples of analytically solvable Hamiltonians are given. 

In the current paper, we present some new examples. Since we are studying Hamiltonians of 
the form ([5]), we are dealing with tridiagonal interaction matrices M with constant entries on the 
diagonal. For tridiagonal matrices, an explicit spectral decomposition (|10p can be found by relating 
these matrices to Jacobi matrices of discrete orthogonal polynomials. So it is natural to look for 
new examples in that area. We shall first describe the example of "Krawtchouk interaction", and 
then indicate how to find other examples. 

3 Krawtchouk interaction 

In this section, let us first collect some (known) properties of Krawtchouk polynomials, and then 
use these to describe the spectrum of a Hamiltonian with a Krawtchouk interaction term. For a 
list of hyper geometric orthogonal polynomials, see [25] or [26]. 

3.1 Krawtchouk polynomials 

For a fixed positive integer parameter and a real parameter p (0 < p < 1), the Krawtchouk 
polynomial of degree i (i = 0, 1, . . . , A^) in the variable x is defined by [25-27] 

Ki(x) ^ K,{x;p,N) = 2F1 (";^;^';^) • (28) 

Herein, 2F1 is the usual Gauss hyper geometric series [28] 

In (j28p . the series is terminating because one of the numerator parameters is a negative integer. 
Note that in (|29p we use the notation of the raising factorial, which can also be rewritten by means 
of a (generalized) binomial coefficient: 

(a)fc = a{a + l)---{a + k-l) = {-if {^^^^ k\ 

The Krawtchouk polynomials satisfy a discrete orthogonality relation of the form 

N 

'^w{x)Ki{x)Kj{x) = hi5ij, (30) 

x=0 

where w{x) is a weight function in x and hi is a function depending on i: 

wix) = (^^^^p^l-pf-^ (x = 0,l,...,iV); /i, = _l_(^i_Zy. (31) 
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Recall that the recurrence relation for Krawtchouk polynomials is given by 

- xKiix) = i(l - p) Ki_iix) - [p{N i(l - p)] Ki{x) + p{N - i) Ki+^{x). (32) 

For future purposes we will however be interested in an orthonormality condition, so we define the 
orthonormal Krawtchouk polynomials by 



Ki{x) = k,{x-p,N) 



yjw{x) Ki{x) 



i = 0,l,2, 



(33) 



Now we can state the following property: 

Lemma 1 Let Mk he the tridiagonal (N + 1) x (N + l)-matrix 

( Fo -El 

-El El —E2 

Mk 







V 



-E2 F2 








—En 
-En En ) 



(34) 



where 



Ei = ^/p{l-p)yJi{N -i + l), Ei = Np+{l- 2p)i, 
and let U be the {N + 1) x (N + l)-matrix with matrix elements 



U,, = Ki{j) = 
where i, j = 0,1, . . . , N . Then 
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p'+^{l-p) 



N-i-j 



1/2 min(ij) 

E \k)\k) ( '^^^k 

^ p' ' 

A;=0 \kJ ^ 



UU^ = U^U = I 



and 



Mk = UDU^ 



(35) 



(36) 



(37) 



where D = diag{0, 1,2..., N). 
Proof. We have that 



N N 

)ij ^ X] ^ikUjk = ^ Ki{k)Kj{k) = 6ij 

k=0 k=0 

by the orthogonality relations (I30p . So UU'^ = I, hence U'^ is the inverse of U and thus also 
U'^U = /. Furthermore, notice that 



E^ = i{l-p) 



hi 



and Ei^i = p{N — i) 



i+l 



hi 



Equation (j32p can then be rewritten as a recurrence relation for the orthonormal Krawtchouk 
polynomials Ki{x): 

xKi{x) = -Ei Ki^i{x) + Ei ki{x) - E,+i ki+i{x). (38) 
Then we have, using Eq = En+i = 0, 



N 



{MKU)ij = ^{MK)ikUkj 



k=0 



-Ei ki_i{j) + E, ki{j) - Ei+i k,+i{j) 

j k,{j) = {UD)ij, 
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so MkU = UD or Mk = UDU^. □ 
So we now have a good candidate interaction matrix Mk- In order to describe systems of the 
form ([5]), however, the diagonal entries Fi of Mx should be constants (i.e. independent of i). We see 
from (I35p that this is the case for p = 1/2. So this leads to a new analytically solvable Hamiltonian 
of the form ([5]). 

3.2 Hamiltonian with Krawtchouk interaction 

Consider a linear chain of n identical harmonic oscillators, with a nearest neighbour interaction 
that is given by 

n o n— 1 

\ - / pi. muj cm \ - ^- „ „ , , 

= E + ^^') - 1" E VK^-O 'Ir'ir+l- (39) 

r=l r=l 

We shall refer to the interaction term as "Krawtchouk interaction". The purpose is to find the 
analytic solution for the spectrum of Hk- It is easy to see that this Hamiltonian can be written in 
matrix form, like (|9l): 



= — ( p\ ■■■ pi 
2m \ 



( Pi 



\ Pn 



( Qi 



\ (In 



(40) 



where Mk is the matrix ()34p with = n — 1 and p = 1/2. But for the matrix Mk we have 
an explicit spectral decomposition, given in Lemma [H and the eigenvalues of Mk are given by 
0, 1, . . . , n — 1. Using this result, and following the general procedure described in section [21 one 
introduces here the following quantities: 



^2 _ c(n_jO ^ _ 1) ^ _ _ 2j + 1), (j = 1, 2, . . . , n). (41) 

The interaction matrix (tj^ — '^^"^ '*''' )/ + cMk is positive definite if all quantities under the square 

root symbol are positive. Since c (and w^) is positive, lo"^ — ^^"2 + c(i — 1) (j = 1,2,..., n) is 
an increasing sequence as j increases. So this condition leads to c < 2w^/(n — 1), or the "coupling 
strength" should be sufficiently small. Now we have: 

Proposition 2 The Hamiltonian Hk given by (I39p is analytically solvable. The explicit spectrum 
of Hk follows from ([Ml) ■ 

" 1 

HK\ki, ...,kn)= '^hxvjikj + -) \ki, . . .,kn), (42) 



where the constants tOj are given by ujj = \J iJ^ — c{n — 2j + l)/2. 

Finally, notice that the interaction term in (j39p is invariant under the reflection r ^ n — r. 



3.3 Remark 

It is clear that the general procedure worked out here for the Krawtchouk polynomials works in 
general for discrete orthogonal polynomials. So in order to find other interesting examples, one 
can go through the list of discrete orthogonal polynomials [25-27] and their (/-analogues. The basic 
restriction, in order to have Hamiltonians of the form ([5]), is that the diagonal elements in the 
interaction matrix are constant (for specific values of the parameters). An investigation of this 
restriction has shown that, apart from the Krawtchouk polynomials with p = 1/2, there are only 
the following cases to be considered: the Hahn polynomials with (3 = a and the dual g-Krawtchouk 
polynomials with c = —1. We shall now study these cases and the corresponding Hamiltonians. 
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4 Hahn interaction 



4.1 Hahn polynomials 



Hahn polynomials are again a class of discrete orthogonal polynomials Qi{x; a, l3, N), characterized 
by a positive integer parameter and two real parameters a and (3 (for orthogonality, one should 
have a > — 1 and /? > —1, or a < —N and /? < —N). The Hahn polynomial of degree i 
{i = 0,1, . . . , N) in the variable x is defined by [25, 26] 



Qiix) = Qi{x;a,f3,N) 



3F2 



+ a + P + 1, —X 
a + l,-N 



(43) 



where 3-F2 is the generalized hypergeometric series (which is terminating here due to the numerator 
parameter —i). The Hahn polynomials satisfy a discrete orthogonality relation: 



N 



'^w{x)Qi{x)Qj{x) = hidij. 



x=0 



where 



w{x) 
hi 



(x = 0,l,...,iV); 



a + x\l N + l3 - x'^ 
x J \ N — X 
i\{N -i)\{i + a + f3 + l)iv+i(/3 + 1) 



iV!2 {2i + a + 13 + l){a + l)i 
The recurrence relation for Hahn polynomials is given by 

- xQi{x) = Ai Qi+i{x) - {Ai + d) Qiix) + d Qi-i{x), 

where 

+ a + + 1)(^ + Q + 1)(A^ - i) i{i + a + (3 + N + + (5) 

"~ {2i + a + 15 + l){2i + a + (3 + 2) ' ' 

The orthonormal Hahn polynomials are defined by 

^Jw{x) Qi{x) 



{2i + a + (5){2i + a + (5 + 1) 



Qi{x) 



0,1,2,... ,N. 



Then we have a similar lemma as Lemma [TJ 
Lemma 3 Let Mq be the tridiagonal (N + 1) x + l)-matrix 

( Fo -El 
-El Fi —E2 



\ 







-Eo 







\ 







—En 
-En Fn J 



where 



Ei 



H{i + a){i + f5){i + a + f5){i + a + f5 + N + 1){N -i + 1) 
{2i + a + f5)'^{2i + a + (5 - l){2i + a + 15 + 1) 



^ _ iV (g - 0) [{a + f5){N - 2i) - 2i{i + 1)] 
* ~ Y 2(2i + a + f5){2i + a + j5 + 2) 



(44) 



(45) 



(46) 



(47) 



(48) 
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and let U be the {N + 1) x {N + l)-matrix with matrix elements 

Uij = Qi{j) (49) 

where i, j = 0,1, . . . , N . Then 

UU'^ = U'^U = I and Mq = UDU'^ (50) 
where D = diag{0, 1,2... ,N). 

Proof. The proof is essentially the same as that of Lemma [1] 

N N 

{UU^)ij = UikUjk = Yl Q^{k)QJ{k) = 

k=0 k=0 

by the orthogonality relations (j44p . hence UU'^ = U'^U = I. Furthermore, equation (j45p can then 
be rewritten as a recurrence relation for the orthonormal Hahn polynomials Qi{x): 

xQi{x) = -Ei Qi-i{x) + Fi Qi{x) - Sj+i QiJ^i{x), (51) 

and this implies MqU = UD or Mq = UDU^. □ 
It remains to be determined when the diagonal of (j47p is constant, in other words when Fi is 
independent of i. Following (148p . this happens when (3 = a. Note that in that case the parameter 
should satisfy q > — 1 or a < —N. It is worthwhile mentioning that there are two other cases 
where the diagonal of (j47p is "almost constant" : 

• When j3 = -a (with -1 < a < 1) one finds that all Fi = {N - a)/2 for i 
Fo = N{a + l)/2. 

• When P = -2N -2 -a (with -2 - N < a < -N) one finds that all Fi 
i = 0,1,..., N -1, hut Fn = N{N + a + 2)/2. 

For these cases, it could still be interesting to consider the Hamiltonian ([9]) built from the corre- 
sponding interaction matrix. The Hamiltonian, however, is not of the form ^ as either the first 
or last oscillator in the chain would play a special role and give rise to an extra term (either in qf 
or else in q'^); so the chain would no longer consist of identical oscillators, but have one of them 
different from the others. 

The Hamiltonians corresponding to P = a will now be considered in the next subsection. 



= 1,2, ... ,A^, but 
= -(a + l)/2 for 



4.2 Hamiltonian with Hahn interaction 



Consider a linear chain of n identical harmonic oscillators with a nearest neighbour interaction 
given by 



He 



V2m 



muj 



cm 
~2 



E 



/ r(n — r)(r + 2a)(r + 2a + n) 
(2r + 2a-l)(2r + 2a + l) 



Irqr+l, 



(52) 

_ _ II I zr' -I- zfv — ] ]i -t- zry -i- i i 

r=l r=l 

the interaction term to be referred to as the "Hahn interaction", and where a is some parameter 
with a > — 1 or a < — n + 1 (guaranteeing that the expression under the square root is positive). 
This Hamiltonian can be written in matrix form: 



He 



1 
2m 



( Pi 



\ Pr. 



+ - 



m 



Qi 



where Mq is the matrix ([^7|) with P = a and A'" 
spectral decomposition of Mq is again diag(0, 1,2. 



qi 



n 
,n 



(C 



c{n — 1) 



)I+cMq) 



( ft 



\ (In 



(53) 



1. Since the diagonal matrix D in the 
1), it follows that: 
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Proposition 4 The Hamiltonian Hq given by (|52p is analytically solvable. The spectrum of Hq 
is exactly the same as that of Hk, and given by (j42p and (|4ip . 



Also the condition for positive definiteness is the same as in the Krawtchouk case, namely c < 
luP' jin — 1). Note that for q = 1/2 the form of the interaction term is considerably simpler: 



ry(i) \^ f Pr , m-^^ ~2\ cm ^ \/(n-r)(n + r + l) „ ^ 



V2m 2 "'^y 2 

r=l r=l 



Observe that for a — > +oo, the interaction (|52p reduces to the Krawtchouk interaction (|39p . Also 
note that under the reflection r — > n — r in the interaction term, the Hamiltonian reads 

-Z^l^ + ^^rj-^Z. ^ ^r^r+l. (55) 

r=l r=l 

This can also be obtained by taking a = —n — 1/2 in (j52]) . 



5 g-Krawtchouk interaction 

5.1 The dual q'-Krawtchouk polynomials 

For a fixed positive integer parameter A^, and real parameters g > and c < [§, the dual q- 
Krawtchouk polynomial of degree i in the variable X{x) = + cq^~^ is defined by [25, 26] 

' q-N Q 1 , (56) 

where 3(^)2 is the g-generalized (or basic) hypergeometric series (which is terminating here due to 
the numerator parameter q~^). Recall that [28] 

f a,b,c \ ^ {a;q)k{b;q)k{c;q)k k 
392 , ; 9; ^ — > , n — — — > 



' ' / ^Qid;q)k{e;q)k{q;q)k 



where the g-Pochhammer symbol is [28] 

(a; q)n = (1 - a)(l - ag) • • • (1 - ag"-^). 
The dual g-Krawtchouk polynomials satisfy the discrete orthogonality relation: 

N 

J2Hx)K^{X{xy,q)K,{X{xy,q) = h^j, (57) 



x=0 



where 



^(x) = : ;g;;~^),-(-.)g.(2iv-.) (, = o, i, . . . , at); 



(g;g)^(cg;g)^(l - cq 



standard literature, this parameter is usually denoted by c, but we replace it by c in order not to confuse with 
the notation for the coupling constant c. 
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The recurrence relation for dual (/-Krawtcliouk polynomials is given by 

- (1 - q--){l - cq--^)K,{X{xy,q) = cq'^ {1 - q')K;,^{\{x);q) 

- [(1 - g-^) + -cq-^{l - q')]K,{\{x)-, (?) + (!- g*-^)K.+i(A(x); q) 

The orthonormal dual g-Krawtchouk polynomials are defined by 



Ki{\{x)-q) 
Then we have the usual lemma: 



yJw{x)Ki{\{x)]q) 



//i7 



i = 0,l,2, 



(58) 



(59) 



Lemma 5 Let M^q he the tridiagonal {N + 1) x (N + l)-matrix 

( Fo -El 
—El Fi —E2 
MKq= -E2 F2 ■■• 



V 







—En 
-En En J 



(60) 



where 



E, 



cq-^{l - q'){l - g^-l-^) = (1 - q'-^) + cq-^ {1 



and let U be the {N + 1) x {N + l)-matrix with matrix elements 

Uij=k,{\{j)-q) 

where i,j = 0,1, . . . , N . Then 

UU'^ = U'^U = I and Mxq = UDU^ 
where D = diag{{l - q-^){l - cq^'^)), {j = 0,1,2 ... ,N). 

The proof is the same as those given before, and uses the orthonormality of Ki{\{j); q) and 
(1 - (7-")(l - cg"-^)Ki(A(x); q) = -Ei Ki.i{X{x); q) + F, ki{\{x);q) - Ei+i Ki+i(A(x); q). (64) 



(61) 
(62) 
(63) 



This case is interesting because for c = — 1 the diagonal of (j60p is constant, in other words then 
Fi is independent of i. This leads again to a Hamiltonian of the type i^. 

5.2 Hamiltonian with dual g-Krawtchouk interaction 

Now we consider a linear chain of n identical harmonic oscillators with a nearest neighbour inter- 
action given by 



r=l r=l 



<ir)-^Y. 2Vg'^+^-2"(l-(?-)(l-g"--) qrqr+l, 



n-1 



(65) 



the interaction term to be referred to as "dual g-Krawtchouk interaction", where q > 0. This 
Hamiltonian can be written in matrix form: 



1 

2m 



^Kq = — \P\ ■■■ Pn 




( Qi 



q[ ■■■ 



qi {{uj^-c{l-q'-^))I + cMKq) 



\ 



(66) 
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where Mj^g is the matrix (|60p with c = — 1 and N = n — 1. It fohows that the spectrum of is 
given by (|12]) . with 

= _ c(i _ gi-n) + c(l - + = -^^2 + c(gJ-"+i - g-J), (j = 1,2,... 

(67) 

For positive definiteness of the interaction matrix, all quantities under the square root must be 
positive. It is easy to see that uj'^ + c(g-'~""*"^ — q~^) (j = 1,2, ... ,n) is an increasing sequence 
of j when q > 1 and a decreasing sequence of j when q < 1. So for q > 1 the condition means 
c < gu;2/(l — g^""), while for < g < 1 we need c < g"'u;2/(l — 5"-+^). To conclude, we have 

Proposition 6 The Hamiltonian H^q given by (j65p is analytically solvable. The explicit spectrum 
of Hkq is given by (p^ , where the constants ujj are given by ujj = \J u"^ ^ c{q^^"'^^ — q^^). 

6 Some properties of the spectra and conclusion 

The spectrum of each of the Hamiltonians given here is of the form (j24p . thus it is discrete but 
infinite dimensional. In order to appreciate the differences of the various examples given here, we 
shall plot the energy levels of the singly excited states (the single phonons, or the simple vibrations) 
of the system. These are the levels of the n states |1, 0, . . . , 0), |0, 1, 0, . . . , 0), . . . , |0, . . . , 0, 1) (in 
the notation of (|24|) ). So, following (j24|) . these levels are given by 

Eq + huJi,Eo + hu!2, ■ ■ ■ ,Eo + tuvn, (68) 

where 

1 " 

Eo = -Y,f^j- (69) 
i=i 

In order to illustrate the spacing of the energy levels of the singly excited states, it is sufficient to 
plot the values of (a;i,a;2, • • • ,Wn)- We plot these values in Figure 1, for n = 12, in four different 
cases: 

(a) The Hamiltonian (jl]) with constant nearest neighbour interaction, where the values iVj are 
given by ([27|) . plotted in Figure 1(a). 

(b) The Hamiltonian with Krawtchouk interaction (j39p or with Hahn interaction (j52p . which have 
the same spectrum and where the values uij are given by (f4T|) . plotted in Figure 1(b). 

(c) The Hamiltonian with (7-Krawtchouk interaction (jGSp where q > 1, for which the ujj are given 
by (i67|) . plotted in Figure 1(c). 

(d) The same Hamiltonian (j65p but with q < 1, for which the ujj are also given by (j67l) . plotted 
in Figure 1(d). 

The values of c and q are appropriately chosen (see the figure caption for actual values) in order 
to illustrate the typical energy level spacing properties for each case. 

For a Hamiltonian with constant nearest neighbour interaction like ([2]) or ([3]), the levels are 
wider apart in the middle of the spectrum, and closer to each other near the top and the bottom 
of the spectrum, see Figure 1(a). The property is known, and was e.g. also observed in [21]. For 
a Hamiltonian with a Krawtchouk interaction or a Hahn interaction like (j39p or (j52p . the energy 
level spacing decreases as the energy increases, a phenomenon also typical for molecular spectra, 
see Figure 1(b). Finally, for a Hamiltonian with a g-Krawtchouk interaction like (|65p . the energy 
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level properties depend on whether 0<g<lorg>l. For q > 1, one observes just the opposite 
of a constant nearest neighbour interaction: the energy level spacing is small near the middle of 
the spectrum, and larger near the top and bottom of the spectrum (but they are wider apart near 
the top than near the bottom), see Figure 1(c). For < g < 1, the levels behave similarly: the 
energy level spacing is small near the middle of the spectrum, and larger near the top and bottom 
of the spectrum, but now they are wider apart near the bottom than near the top, see Figure 1(d). 
It should be noted that for the cases (a), (b) and (c) the order of the levels from bottom to top 
correspond to the states |1, 0, . . . , 0), |0, 1, . . . , 0), . . . , |0, 0, . . . , 1) in this order, whereas for (d) it is 
just the opposite. This is related to the fact that the sequence ujj {j = 1,2, . . . , n) is an increasing 
sequence in the cases (a), (b) and (c), but a decreasing sequence in the case (d) (see the remark 
following eq. ([67|) ). 

To conclude, we have presented a number of Hamiltonians for a quantum system consisting of 
a linear chain of identical oscillators with a nearest neighbour interaction term depending on the 
position in the chain, in general given by ([5]). Although such Hamiltonians are always numerically 
solvable, we have focused on the question when such a Hamiltonian is analytically solvable. The 
connection between the interaction matrix of the matrix expression of the Hamiltonian and Jacobi 
matrices of discrete orthogonal polynomials has led to new analytically solvable Hamiltonians, with 
interesting spectral properties. 

Our original interest in Hamiltonians with a nearest neighbour interaction of the form ([2]) 
stems from the fact that this Hamiltonian can also be quantized in an alternative way, namely 
as a Wigner quantum system [29-31], leading in particular to finite spectra and non-commutative 
coordinates [19,20,32-35]. Our aim is to investigate when Hamiltonians with a general interaction 
matrix, of the form ([9]), can still be solved as a Wigner quantum system, and investigate the 
properties of such a solution. We hope to report on that topic in the near future. 
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Figure 1: Energy levels of the n single phonon states |0, . . . , 0, 1, 0, . . . , 0), for n = 12. In each of 
the cases, h = lo = 1. The four cases correspond to: (a) a Hamiltonian ^ with constant nearest 
neighbour interaction [for c = 0.5]; (b) a Hamiltonian with Krawtchouk interaction (j39p or with 
Hahn interaction (I52p [for c = 0.18]; (c) a Hamiltonian with g-Krawtchouk interaction ()65p where 
g > 1 [here, q = 1.6 and c = 1.0]; (d) the same g-Krawtchouk interaction (j65p where < g < 1 
[here, q = 0.7 and c = 0.01]. The levels are rescaled, so that the lowest and highest levels match in 
the four cases. 

(a) (b) (c) (d) 
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